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We characterize the thermodynamical equihbrium states of axisymmetric Euler-Beltrami flows. 
They have the form of coherent structures presenting one or several cells. We find the relevant 
control parameters and derive the corresponding equations of state. We prove the coexistence of 
several equilibrium states for a given value of the control parameter like in 2D turbulence [Chavanis 
& Sommeria, J. Fluid Mech. 314, 267 (1996)]. We explore the stability of these equilibrium states 
and show that all states are saddle points of entropy and can, in principle, be destabilized by a 
perturbation with a larger wavenumber, resulting in a structure at the smallest available scale. This 
mechanism is therefore reminiscent of the 3D Richardson energy cascade towards smaller and smaller 
scales. Therefore, our system is truly intermediate between 2D turbulence (coherent structures) and 
3D turbulence (energy cascade). We further explore numerically the robustness of the equilibrium 
states with respect to random perturbations using a relaxation algorithm in both canonical and 
microcanonical ensembles. We show that saddle points of entropy can be very robust and therefore 
play a role in the dynamics. We evidence differences in the robustness of the solutions in the 
canonical and microcanonical ensembles. A scenario of bifurcation between two different equilibria 
(with one or two cells) is proposed and discussed in connection with a recent observation of a 
turbulent bifurcation in a von Karman experiment [Ravelet et aL, Phys. Rev. Lett. 93, 164501 
(2004)]. 

PACS numbers: 



I. INTRODUCTION 

The statistical mechanics of systems with long-range interactions has recently attracted a lot of attention [l| . Typical 
systems with long-range interactions include self-gravitating systems [H, Q , two-dimensional vortices 0, non- neutral 
plasmas [H , free electrons lasers [6[ and toy models such as the Hamiltonian Mean Field (HMF) model [3, [H . Unusual 
properties of systems with long-range interaction such as negative specific heats or ensembles inequivalence have been 
evidenced and Hnked with lack of additivity [9|]. In addition, a striking property of these systems is the rapid formation 
of quasi stationary self-organized states (coherent structures) such as galaxies in the univers e fiol, l arge scale vortices 
in geophysical and astrophysical flows |lll-[l5| or quasi-stationary states in the HMF model [l6l-[l8l|. These QSSs can 
be explained in terms of statistical mechanics using the theory developed by Lynden-Bell [T9| for the Vlasov equation 
or by Miller [20| and Robert & Sommeria [21] for the 2D Euler equation. 

Two-dimensional vortices interact via a logarithmic potential. Interaction of vortices in 3D turbulence is weaker 
than in 2D turbulence, but still long-range. Due to dissipative anomaly and vortex stretching, statistical mechanics 
of 3D turbulence has so far eluded theories. Recent progress was recently made considering 3D inviscid axisymmetric 
flows [22I, I23I that are intermediate between 2D and 3D flows: they are subject to vortex stretching like in 3D 
turbulence, but locally conserve a scalar quantity in the ideal limit, like in 2D turbulence. It is therefore interesting 
to study whether these systems obey the peculiarities observed in other systems with long-range interactions such as 
violent relaxation, existence of long-lived quasi-stationary states, negative specific heats and ensembles inequivalence. 

The general study of the stability of axisymmetric flows, and the possible occurrence of phase transitions, is difficult 
due to the presence of an infinite number of Casimir invariants linked with the axisymmetry of the flow. In a previous 
paper [2^, hereafter Paper I, we have considered a simplified axisymmetric Euler system characterized by only three 
conserved quantities: the fine-grained energy E^-^- ^ the helicity H and the angular momentum /. We have developed 
the corresponding statistical mechanics and shown that equilibrium states of this system have the form of Beltrami 
mean flows on which are superimposed Gaussian fluctuations. We have shown that the maximization of entropy S 
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at fixed helicity iJ, angular momentum / and microscopic energy E^-^- (microcanonical ensemble) is equivalent to 
the maximization of free energy J = S — f3E^-^- at fixed helicity H and angular momentum / (canonical ensemble). 
These variational principles are also equivalent to the minimization of macroscopic energy E^-^- at fixed helicity H 
and angular momentum /. This provides a justification of the minimum energy principle (selective decay) from 
statistical mechanics. We have furthermore discussed the analogy with the simplified thermodynamical approach of 
2D turbulence developed in based on only three conserved quantities: the fine-grained enstrophy T^'^'^ the energy 
E and the circulation F. We have shown that equilibrium states of this system have the form of Beltrami mean flows 
(linear vorticity-stream function relationship) on which are superimposed Gaussian fluctuations. We have shown 
that the maximization of entropy S at fixed energy E, circulation F and microscopic enstrophy T^'^' (microcanonical 
ensemble) is equivalent to the maximization of grand potential S = S — a2T2^' at fixed energy E and circulation F 
(grand microcanonical ensemble). These variational principles are also equivalent to the minimization of macroscopic 
enstrophy F2^" at fixed energy E and circulation F. This provides a justification of a minimum enstrophy principle 
(selective decay) from statistical mechanics. In the analogy between 2D turbulence and 3D axisymmetric turbulence, 
the energy plays the role of the enstrophy. 

In the present paper, we study more closely the equilibrium states of axisymmetric flows and explore their stability. 
We show that all critical points of macroscopic energy at fixed helicity and angular momentum are saddle points, so 
that they are unstable in a strict sense. Indeed, there is no minimum (macroscopic) energy state at fixed helicity 
and angular momentum (either globally or locally) because we can always decrease the energy by considering a 
perturbation at smaller scales. This is reminiscent of the Richardson energy cascade in 3D turbulence. Inversely, 
in 2D turbulence, there exists minimum enstrophy states that develop at large scales (inverse cascade). Therefore, 
our system is intermediate between 2D and 3D turbulence: there exists equilibrium states in the form of coherent 
structures (that are solutions of a mean field differential equation) like in 2D turbulence, but they are saddle points of 
macroscopic energy and are expected to cascade towards smaller and smaller scales like in 3D turbulence. However, 
we give arguments showing that saddle points can be robust in practice and play a role in the dynamics. Indeed, they 
are unstable only for some particular (optimal) perturbations and can persist for a long time if the system does not 
spontaneously generate these perturbations. Therefore, these large-scale coherent structures can play a role in the 
dynamics and they have indeed been observed in experiments of von Karman flows [25]. In order to make this idea 
more precise, we have explored their stability numerically using phenomenological relaxation equations derived in 
[23| . We have found some domains of robustness in the parameter space. In particular, the one cell structure is highly 
robust for large values of the angular momentum I > Ic and becomes weakly robust for low values of the angular 
momentum. In that case, we expect a phase transition (bifurcation) from the one-cell structure to the two-cells 
structure. We have also found that the value of the critical angular momentum Ic changes depending whether we 
use relaxation equations associated with a canonical (fixed temperature) or microcanonical (fixed microscopic energy) 
description. At low temperatures T, we have evidence a new kind of "ensembles inequivalence" characterizing the 
robustness of saddle points with respect to random perturbations. 

The paper is organized as follows: In Sec. [Ill we set-up the various notations and hypotheses we are going to 
use. The computation and characterization of equilibrium states is done in Sec. Illli The stability analysis of these 
equilibrium states is performed in Sec. IIII Gl where we show analytically that all states are unstable with respect to 
large wavenumber perturbations. We evidence a process of energy condensation at small scales that is reminiscent of 
the Richardson cascade. We explore numerically the robustness of the equilibria in both canonical and microcanonical 
ensembles in Sec. llVi Our numerical method is probabilistic and rather involved. A discussion of our results is done 
in Sec. [Vl where a bifurcation scenario relevant to the turbulent experimental von Karman flow is suggested. 



II. THEORETICAL SET-UP 
A. Hypotheses and Notations 



We consider a system with a cylindrical geometry enclosed in the volume delimited above and below by surfaces 
z = and z = 2/i, and radially hy < r < R. Like in Paper I, we consider an axisymmetric Euler-Beltrami system 
characterized by a velocity field u, with axisymmetric time averaged u. We furthermore assume that the only relevant 
invariants of the axisymmetric Euler equations for our problem are the averaged energy E = ^ J dr, the averaged 
hehcity H = J u ujdr and the averaged angular momentum / = Jadr where a = rue- We introduce the potential 
vorticity = uoq/v and the stream function such that Ur = —r~'^dzip and Uz = r~'^dril^. They are related to each 
other by the generalized Laplacian operator 

A = — ^ + i— (-^\ - (1) 

* dz"^ r dr \r dr J 
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In actual turbulent von Karman experiments, we have been able to observe that the largest part of the kinetic energy is 
contained in the toroidal motions. It is therefore natural, as a first elementary step, to consider a model in which only 
toroidal fluctuations are considered, and suppose that the fluctuations in the other (poloidal) directions are simply 
frozen. With such an assumption, poloidal vorticity fluctuations are allowed, but toroidal vorticity fluctuations are 
excluded. We therefore only include a fraction of the vorticity fluctuations, that presumably become predominant at 
small scale, due to the existence of vortex stretching. As shown below and in the next paper [30|, this simplification 
however still allows for vortex stretching and energy cascades towards smaller scales, and leads to predictions that are 
in good agreement with experiments. Moreover, our hypotheses lead to a model that is self-contained and analytically 
tractable. According to our hypotheses, neither ^ nor ^jj fluctuates in time: i = i and ?/; = ?/;. In that case, the 
conserved quantities can be rewritten 

^'^■-\m^\(^), (2) 

H={la), (3) 
I={a), (4) 



where (/) denotes the spatial average [32| 
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The helicity and the angular momentum are robust constraints because they can be expressed in terms of coarse- 
grained quantities ^ and a. By contrast, the energy is a fragile constraint because it cannot be expressed in terms of 
coarse-grained quantities. Indeed, it involves the fluctuations of angular momentum a^. To emphasize that point, we 
have introduced the notation E = E^-^- to designate the fine-grained (microscopic) energy. Splitting a into a mean 
part a and a fluctuating part Sa^ we define the coarse-grained (macroscopic) energy by 

Then, the energy contained in the fluctuations is simply 

Ef,^^,^Ef-^--E-^- = \{^), (7) 

where 

a2=^-a^ (8) 

is the local centered variance of angular momentum. We stress that the microscopic energy E = E^-^- is conserved 
while the macroscopic energy E^-^- is not conserved and is likely to decrease (see below). 

B. The two statistical ensembles and the selective decay principle 

In Paper I, we have developed a simplified thermodynamic approach of axisymmetric flows under the above- 
mentioned hypothesis. Let p(r, r]) denote the PDF of a and let us recall the expression of the entropy 

S[p] = - J plnpdrdr]. (9) 

We have proven the equivalence between the microcanonical ensemble 



max 



^x{S[p]\Ef-^-,H,I, J pdr] = l}, (10) 
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and the canonical ensemble 

max{ J[p] = S- pE^-^- \ H, /, / pdr] = 1}. (11) 
P,^ J 

In each ensemble, the critical points are determined by the first order condition SS — /3SE^-^- — jiSH — adi = 0. The 
equilibrium distribution is Gaussian 



1/2 



Kr,r?)=(^) e-^^^-^K (12) 



the mean flow is a Beltrami state 
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V', (13) 

?=-4--, (14) 
and the centered variance of angular momentum is 

^2 = J. (15) 

These equations determine critical points of the variational problems (flO]) and (fTT]) that cancel the first order variations 
of the thermodynamical potential. Clearly, (p!Q|) and (pTj) have the same critical points. Furthermore, it is shown 
in Paper I that (p!Q)) and (pTj) are equivalent for the maximization problem linked with the sign of the second order 
variations of the thermodynamical potential: a critical point determined by Eqs. (p!2]) - (p!5]) is a maximum of S at fixed 
microscopic energy, helicity and angular momentum iff it is a maximum of J at fixed helicity and angular momentum. 
This equivalence is not generic. We always have the implication (pTj) (p!Q|) but the reciprocal may be wrong. Here, 
the microcanonical and canonical ensembles are equivalent due to the quadratic nature of the microscopic energy ~ . 
We note that, according to Eq. (fT5]) . P is positive. In the canonical ensemble, P is prescribed. In the microcanonical 
ensemble, /3 is a Lagrange multiplier that must be related to the energy E = E^-^-. According to Eqs. ([7]) and (fT5]) . 
we find that /3 > is determined by the condition 

Ef-9--E^-s-=Efi^,t = ^. (16) 

This relation shows that T = 1//3 plays the role of a temperature associated with the fluctuations of angular momentum 
(33| . Finally, we have proven in Paper I that the two variational problems (p!Q|) and (pT]) are equivalent to 

max{J[a,?] = -/3^-^-|i7, /}, (17) 



or equivalently 



min{E-^-[a,e]|i7, /}, (18) 



in the sense that the solution of (fTQ|) or (fTTj) is given by Eq. (fT2)) where (a, ^) are the solutions of (pT|) or (fTSj) . 
This justifies a selective decay principle from statistical mechanics. Indeed, it is often argued that an axisymmetric 
turbulent flow should evolve so as to minimize energy at fixed helicity and angular momentum. In general, this 
phenomenological principle is motivated by viscosity or other dissipative processes. In our approach, it is justified by 
the maximum entropy principle (fTOj) of statistical mechanics when a coarse-graining is introduced. In the sequel, we 
shall study the maximization problem (fT8|) since it is simpler than (p!Q|) or ([TTj), albeit equivalent. 

Remark: although the variational problems (fTQ|) and (fTTj) determining equilibrium states are equivalent, this does 
not mean that the relaxation equations associated with these variational problems are equivalent. To take an analogy, 
the Boltzmann (microcanonical) and the Kramers (canonical) equations have the same equilibrium states -the Maxwell 
distribution- but a different dynamics. In the following, we will show that the equilibrium variational problems (fTQ|) 
and (pTj) have no solution. Indeed, there is no maximum of entropy at fixed E^-^-, H and / and no minimum of 
free energy at fixed H and /. All the critical points of (p!Q|) and (pTj) are saddle points of the thermodynamical 
potentials. Then, the idea is to consider the out-of-equilibrium problem, introduce relaxation equations and study the 
robustness of saddle points with respect to random perturbations. For what concerns the out-of-equilibrium problem, 
the microcanonical and canonical ensembles may be inequivalent. We will see that they are indeed inequivalent. 
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III. MINIMUM ENERGY STATES 

A. Critical points 

In this section, we shall study the minimization problem 

mm{E'-^[a,^]\H,I}. (19) 

The critical points of macroscopic energy at fixed helicity and angular momentum are determined by the condition 

SE""-^- + ^6H + aSI = 0, (20) 

where ja (helical potential) and a (chemical potential) are Lagrange multipliers. Introducing the notations B = — 
and C = —a/ the variations on and 5a lead to 

a = Bilj, (21) 

? = 5^ + C, (22) 

which are equivalent to Eqs. (p!3)) - (p!4|) up to a change of notations. In the following, it will be convenient to work 
with the new field (j) = t/j/r. It is easy to check that 

A.^ = i (a0 - ^) , (23) 

where A is the usual Laplacian. Therefore, Eq. ([1]) becomes 

-A0+^=r?, (24) 



and the previous equations can be rewritten 



a = Brcp, (25) 



e = — + C, (26) 

r 



- A(/)+ J =52(/) + Cr, (27) 

with = on the boundary. This is the fundamental differential equation of the problem. Note that a particular 
solution of this differential equation is 



where 6 is solution of 



Vpart = --f^r, (28) 



but it does not satisfy the boundary conditions. Using Eqs. ([3]) and (|4]), the hehcity and the angular momentum are 
given by 



H -CI = B^{(t)^), (29) 



I = B{(l)r). (30) 



These equations are relationships between {B^C) and (H^I). 

Remark: we have not taken into account the conservation of circulation F = (^) because this would lead to a term 
AB/r in the r.h.s. of Eq. (|27|) that diverges as r ^ 0. 
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B. The different solutions 



To construct the different solutions of Eq. ([27]) and study their stability, we shall follow the general procedure 
developed by Chavanis & Sommeria [26| for the 2D Euler equation. We first introduce an eigenmode decomposition 
to compute all critical points of ([19]). Then, we investigate their stability by determining whether they are (local) 
minima of macroscopic energy or saddle points. 



1. The eigenmodes 

We first assume that 

C = 0. (31) 

In that case, the differential equation (|27|) becomes 

-A^^^=B^, (32) 

with ^ = on the domain boundary. We introduce the eigenfunctions of the operator £ = — A + They are 
defined by 

>C0mn = -A0^n + = 5^^0mn, (33) 

with (j)rnn = ou the domain boundary. It is easy to show that the eigenvalues Amn of C are positive (hence the 

notation Amn = ^mn)- Indeed, we have {(l)mn^(t)mn) = Amn((/>mn) and {(l)mnC(t)mn) = ((V(/)mn)^) + (^p) > 0, which 
proves the result. It is also easy to show that the eigenfunctions are orthogonal with respect to the scalar product 

-| nR n2h 

(f9) = j^J^ rdrdz f{r,z)g{r,z). (34) 

Finally, we normalize them so that {(l)mn(l)m' n') = ^mm'^nn'- 

The eigenvalues and eigenfunctions of the operator C can be determined analytically. The differential equation ([33]) 
can be rewritten 

3^4, id4> <p 2 

We look for solutions in the form 0(r, z) = f(r)g{z). This yields 

f r f g 

where the sign of the constant has been chosen in order to satisfy the boundary condition = in z = and z = 2h. 
The differential equation for g is readily solved and we obtain 

g{z) = sin(/^2;), (37) 

with 

K = n^, (38) 

where n is a strictly positive integer. On the other hand, the differential equation for / is 

ry" + rf'-r\K^-B')f-f = 0. (39) 
If we define \^ = B'^ — and x = Ar, the foregoing equation can be rewritten 

x^f" + xf + x^f-f = 0. (40) 
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This is a Bessel equation whose solution is 

f = .h{x). (41) 
Now, the boundary condition (l){R) = implies f{XR) = so that 

\R = jim. (42) 
where jim is the m-th zero of Bessel function Ji . In conclusion, the eigenvalues are 



and the eigenfunctions are 



with the normalization constant 



The mode (m, n) corresponds to m cells in the r-direction and n cells in the 2:-direction. We shall distinguish two 
kinds of modes, according to their properties regarding the symmetry IZ with respect to the plane z = h. The odd 
eigenmodes denoted are such that = —(j^'rnn correspond to n even. They have zero mean value in the 

z direction (J^ dz = 0). For example, the mode (1,2) is a two-cells solution in the vertical direction. The even 

eigenmodes denoted (j)'^^ are such that = (j)'^^ and correspond to n odd. They have non zero mean value in 

the vertical direction (L (/)^^ dz ^ 0). In particular, the mode (1,1) is a one-cell solution. 

Returning to Eq. (|3^. this differential equation has solutions only for quantized values of B = Bmn (eigenvalues) 
and the corresponding solutions (eigenfunctions) are 



(46) 



where we have used the helicity constraint ([29|) to determine the normalization constant. Note that Eq. (|29|) implies 
that Bmn and H have the same sign, so that the square root is always defined. Substituting this result in Eq. (|3Q]) . 
and introducing the control parameter 
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A=^, (47) 
we find that these solutions exist only for A = A^^ with 

Amn = ^^(0mnr)^ (48) 

For the odd eigenmodes we have A = and for the even eigenmodes we have A^^ = -^j— {4>mn'^)'^ • 



2. The continuum 



We now assume that C 7^ and define 



6b = g. (49) 



In that case, the fundamental differential equation ()27p becomes 



A</.B + ^ = S20B + r, (50) 



8 



with (j)B = on the domain boundary. We also assume that B ^ Bmn- In that case, Eq. (|5Q]) admits a unique 
solution that can be obtained by expanding (j)B on the eigenmodes. Using the identity 

/ = ^{f(l)mn)(l)mn^ (51) 
mn 

we get 

^-=E^#f^C„. (52) 

mn 'mn ^ 

Of course, can also be obtained by solving the differential equation (|27|) numerically. Note that this solution is 
even since only the even modes are "excited". Substituting Eq. (|^9|) in Eq. (|3Q|). we obtain 

Then, substituting Eqs. (0^]) and dSSJl in Eq. we get 



A= .. (54) 



This equation gives a relationship between B and A. Then, C is determined by Eq. (|53]) . These equations can 
therefore be viewed as the equations of state of the system. They determine the branch formed by the solutions of 
the continuum. Using 

(^-r) = j:^¥^.^ (55) 

mn ^mn ^ 



and 



we obtain 



^^) = E7^#f^' (56) 

mn \^mn ^ ) 



{<t>Br) + B^^l) = ^ ; 7V "2/2 > 0- (57) 

mn y-^mn ^ ) 

This implies that B is of the same sign as A, hence H. Furthermore, A is an odd function of B. In the sequel, we 
shall consider only cases with > 0, i.e. B > and A > for illustration and figures. 
Note that Eq. (|54|) involves the important function 

F{B)^{<PBr) = J2^ff^^. (58) 

mn 'mn 

For A = 0, the inverse helical potential is 5 = or 5 = B^^ where B^^ is any zero of F, i.e. 

F(B(")) = (</.^(„)r)=0. (59) 

For simplicity, we shall call B^ = Bi^^ the first zero of F. This first zero is always between the first and the second 
even eigenmodes (see Appendix A). Its location with respect to the first odd eigenmode B[2 depends on the aspect 
ratio of the cylinder: for h/R > 0.53, we have B^ > B'12 (case L-for Large aspect ratio) while for h/R < 0.53, 
B^ < B'12 (case S-for Small aspect ratio). 
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3. The mixed solutions 



We now consider the case where C 7^ and B = Bmn- For B = B'^^^ we recover the eigenfunction (/)^^ as a hmit 
case. Therefore, the even eigenmodes are hmit points of the main branch. On the other hand, for B = 5^^, the 
solution of Eq. (|5Q]) is not unique. Indeed, we can always add to the solution (|52]) an eigenmode Xmn<\>'mn' This leads 
to the mixed solution 

m'n' ^ 'm'n') \-^mn) 

The "proportion" Xmn of the eigenmode present in the mixed solution is determined by the control parameter A. 
Taking the norm of (j)M and its scalar product with r, we get 



Vrm'n' I 



m'n' I 



Substituting these results in Eqs. (|29|) and (|3Q|) . we find that Xmn is determined by A according to 

A = ^(^^^)' (62) 

with B = Bl^^. These mixed solutions exist in the range < A < A (5^^) and they form a plateau at constant 
B = B'^^. For Xmn +C)0, we recover the odd eigenmode (j)'^^ at A = and for Xmn = 0, the plateau connects the 
branch of continuum solutions. The mixed solutions are therefore symmetry breaking solutions. They can be seen as 
a mixture of a continuum solution and an eigenmode solution, like in situations with different phase coexistence. 



C. The helical potential curve 

In this section, we plot 5 as a function of A. For given / 7^ 0, this curve determines the inverse helical potential 
1/ ji as a function of the inverse helicity 1/H (conjugate variables). It is represented in Figs. [T]and|4]for the cases L 
and S respectively. One sees that, for a given value of the control parameter A, there exists multiple solutions with 
different values of B. We will see in Sec. IIIIDI that, for a given value of A, the macroscopic energy E^-^- decreases 
as B increases. Therefore, low values of B correspond to high energies states and high values of B correspond to low 
energies states. 

Case L: In this case B^ > B[2 and the curve B{K) looks typically like in Fig. [TJ For a given value of A, we have 
different solutions as represented in Fig. [2j The highest energy solution is a one cell solution (continuum branch), 
that we choose to call "direct monopole". The second one is a two vertical cells solution (mixed branch). The cells 
are symmetric for A = but one of the two cells grows for increasing A. The third highest energy solution is another 
one-cell solution (continuum branch) rotating in a direction opposite to that of the highest energy solution. We 
therefore call it a "reversed monopole". We call these three respective branches of solutions "branch 1" and "branch 
2" for the continuum solutions, and "mixed branch" for the mixed solutions. The branches 1 and 2 connect each 
other at A(5('^), the location of the first even eigenmode. A typical sequence of variation of the stream function with 
increasing A on these three branches is given in Fig. [3l One sees that, as we increase A on the mixed branch, the two 
cells solution, with a mixing layer dX z = h continuously transforms itself into a one cell solution, via a continuous 
shift of the mixing layer towards the vertical boundary. 

Case S: In this case B^ < B'12 and the curve B{K) looks typically like in Fig. 21 The highest energy solution is a 
one cell solution (continuum branch), (direct monopole). The second solution is another one-cell solution (continuum 
branch ) rotating in the opposite direction (reversed monopole). The third solution is a two horizontal cells solutions 
(continuum branch). Some stream functions are represented in Fig. [5l 

Remark: there is a maximum value of Kmax above which there is no critical point of energy at fixed helicity and 
angular momentum. In that case, the system is expected to cascade towards smaller and smaller scales since there is 
no possibility to be blocked in a "saddle point" . This is a bit similar to the Antonov instability in stellar dynamics 
due to the absence of critical point of entropy at fixed mass and energy below a critical value of energy ^2|, i3|, i27|, [28| . 
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-0.05 0.05 0.1 0.15 0.2 0.25 0.3 0.35 

A 

FIG. 1: Top: 5 as a function of A for case L (we have taken = 1.4 and h = 1.2). For a given value of A (we have taken 
A = 0.05), the solutions of the continuum are denoted by red circles and the mixed solutions by green circles. The mixed 
solution branches are drawn using dotted lines. One observes multiplicity of solutions: at given A correspond several solutions 
with different B. 




r lA r 1-4 j. 1.4 r 1-4 



FIG. 2: Example of stream function of the four first solutions for A = 0.05. From left to right: B = 0.3767 (direct monopole), 
B = 3.7874 (vertical dipole), B = 4.1633 (reversed monopole) and B = 5.15. Increasing values from blue to red. By convention, 
we call direct (resp. reversed) monopole the one-cell solution with maximal (resp. minimal) inner stream function-see above. 
For simplicity, we show at each point only one solution, corresponding to a given sign of /. The solution corresponding to 
opposite sign of / can be found by a change 0^—0. 

In that case, the system is expected to collapse (gravothermal catastrophe). It is not yet clear whether a similar 
process can be achieved in experiments of turbulent axisymmetric flows. For A > A^ax^ the system could become 
non-axisymmetric ruling out the theoretical analysis. 
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FIG. 3: Stream function (j) along the three branches of solution at A = 0.004, 0.1 and 0.3 from left to right for each branch. 
Top= branch 1, direct monopole; Middle: mixed branch (vertical dipole); Bottom: branch 2, reversed monopole. For simplicity, 
we show at each point only one solution, corresponding to a given sign of /. The solution corresponding to opposite sign of / 
can be found by a change (j) ^ —cj). 

D. The coarse-grained energy 

In the previous section, we have found several solutions with different values of B for each value of the control 
parameter A < Kmax- According to the variational principle (p!9|) , we should select the solution with the minimum 
macroscopic energy. Combining Eqs. (|3|), dH, (|6|), (|2T]) and (|22]) , we obtain the relation 



H-\CI 



(63) 
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FIG. 4: Top: B as a function of A for case S (here R — 1.4 and h — 0.5). The solutions of the continuum are denoted by red 
circles and the mixed solutions bv efreen circles. The mixed solution branches are drawn usinef dotted lines. 
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FIG. 5: Stream function (j) of the four first solutions at A = 0.02 for case S. From left to right: B — 0.6567 (direct monopole), 
B — 5.2033 (reversed monopole), B — 5.8967 and B — 6.8534. For simplicity, we show at each point only one solution, 
corresponding to a given sign of /. The solution corresponding to opposite sign of / can be found by a change (j) —cj). 



For the eigenmodes {C = 0), we find that 



H Br) 



(64) 



Let us consider the odd eigenmodes (j)'^^ that exist for A = only. They are in competition with each other. We 
see that there is no minimum energy state since the energy decreases when (m, n) increase, i.e. when the eigenmodes 
develop smaller and smaller scales. Therefore, the minimum energy state corresponds to the structure concentrated 
at the smallest accessible scale. For the solutions of the continuum, using Eq. (|53|) , the macroscopic energy is 

'^l-^^V (65) 



H B\ 2B{(t)Br) 

We can easily plot it as a function of B (see Figs. [6]and[7|). Combining Figs. [Tl|4l[6]and[7l we see that, for a given value 
of A, the solution with the smallest macroscopic energy corresponds to the highest 5, i.e. to small-scale structures. 
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This is in complete opposition to what happens in 2D turbulence. In that case, the counterpart of the macroscopic 
energy E*^-^- is the macroscopic enstrophy T2^' and the minimum enstroph y st ate corresponds to structures spreading 
at the largest scale. Strikingly, the bifurcation diagram in 2D turbulence [26| is reversed with respect to the present 
one. 

0.5 
0.4 
0.3 
^ 0.2 
0.1 



"012345678 
B 

FIG. 6: Macroscopic energy E^'^ / H as a function of B for case L. The energy of the even eigenmodes are denoted by red circles 
and the energy of the odd eigenmodes by green circles. 

In conclusion, there is no global minimum of macroscopic energy at fixed helicity and angular momentum. We can 
always decrease the macroscopic energy by considering structures at smaller and smaller scales. Since (p^Oj) . (fTTj) and 
(p!8|) are equivalent, we also conclude that there is no global maximum of entropy at fixed microscopic energy, helicity 
and angular momentum. We may note a similar fact in astrophysics. It is well-known that a stellar system has no 
global entropy maximum at fixed mass and energy [5, ^sl, [13, [lH . This is associated to gravitational collapse (called the 
gravothermal catastrophe in the microcanonical ensemble) leading to the formation of binary stars. However, in the 
astrophysical problem, there exists local entropy maxima (metastable states) at fixed mass and energy if the energy is 
sufficiently high (above the Antonov energy). Similarly, we could investigate the existence of metastable states in the 
present problem. However, we will show in Sec. IIII GI that there is no local minimum of macroscopic energy at fixed 
helicity and angular momentum. All the critical points (|2T]) - ([22]) of the variational problem (p!9|) are saddle points! 

0.5 
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FIG. 7: Macroscopic energy E^'^ / H as a function of B for case S. The energy of the even eigenmodes are denoted by red circles 
and the energy of the odd eigenmodes by green circles. 
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E. Chemical potential curve 



In our system, the chemical potential is a = C /B. For given i7, we have to plot a as a function of / (conjugate 
variables). The chemical potential is zero for the eigenmodes. Using the equation of state (|53|) . we can express a for 
the continuum solutions as 



(66) 



where A is expressed as a function of B by Eq. (|54|) . Therefore, Eq. (|66|) gives a/ \/\H\ as a function of B. Ehminating 
B between Eqs. (|66]) and (|54|) . we obtain a/ \/\H\ as a function of A for the continuum. For the mixed solutions, we 
have 



(67) 



corresponding to straight lines as a function of y |A|. The chemical potential curve a/ ^\H\ as a function of 
represented in Fig. [8] for case L and in Fig. [9] for case S. For fixed this gives a as a function of /. 




FIG. 8: Chemical potential versus angular momentum in case L. The chemical potential of the even eigenmodes are denoted 
by red circles and the chemical potential of the odd eigenmodes by green circles. The mixed solution branches are drawn using 
dotted lines. 




FIG. 9: Chemical potential versus angular momentum in case S. The chemical potential of the even eigenmodes are denoted 
by red circles and the chemical potential of the odd eigenmodes by green circles. The mixed solution branches are drawn using 
dotted lines. 
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F. Caloric curve 



If we come back to the initial variational problem ([TOj), the caloric curve should give /3 as a function of the 
microscopic energy E = E^-^- (conjugate variables) for fixed values of H and /. Now, the temperature is determined 
by the expression 

Ef-y-=E'^-^- + ^. (68) 

For given H and /, we can determine the discrete values of Bf^^) and the corresponding discrete values of E"^^^' as 
explained previously. Then, for each discrete value, the temperature is related to the energy by Eq. ([68|) . Therefore, 
the mean flow (Beltrami state) is fully determined by H and / and, for a given mean flow, the variance of the 
fluctuations (temperature) is determined by the energy E = E^'^' according to 

^T = E-El;^y (69) 

In conclusion, the caloric curve T{E)^ or more properly the series of equilibria, is formed by a a discrete number of 
straight lines with value at the origin —E^^^^ and with constant specific heats C = dE/dT = 1/2. The specific heat 
is positive since the microcanonical and canonical ensembles are equivalent in our problem. 



G. Stability analysis 

In this section, we prove that the critical points of macroscopic energy at fixed helicity and angular momentum are 
all saddle points. A critical point of macroscopic energy at fixed helicity and angular momentum is a minimum (resp. 
maximum) iff the second order variations 

5^F = 5^E'-3- - ^ 

= I (^\rS^S<l> + ^ - y^Sa^ rdrdz, (70) 

are definite positive (resp. definite negative) for all perturbations that conserve helicity and angular momentum at 
first order, i.e. SH = {S,Sa) + {o'S^) = and SI = (Sa) = 0. Adapting the procedure of Chavanis & Sommeria [26| to 
the present context, we shall determine sufficient conditions of instability. 

(i) Let us prove that there is no local maximum of macroscopic energy at fixed angular momentum and helicity. 
Consider first the even solutions, including the continuum solutions and the even eigenmodes. We choose a perturba- 
tion such that 6a is odd and 6^, = 6(f) = 0. For symmetry reason, this perturbation does not change / nor H at first 
order. On the other hand, for this perturbation 6'^F = J rdrdz > 0. Consider now the odd eigenmodes. We 
choose a perturbation of the form 6£^ = 6(j) = and 6a = B^r(j)B^ , where (j)B^ is the first continuum solution such that 
(^(/>B^) = 0. For this perturbation, we have {6a) = B^{r(j)B^) = 0, {a6^,) = and {(,6a) = {B'^^)'^B^{(j)B^(t)'mn) ~ ^ 
since (j)'^^ is orthogonal to (pB^- Therefore, this perturbation does not change the helicity and the angular momentum 

at first order. On the other hand, for this perturbation 6'^F = J rdrdz > 0. As a result, the critical points 
of macroscopic energy at fixed helicity and angular momentum cannot be energy maxima since we can always find 
particular perturbations that increase the energy while conserving the constraints. 

(ii) Let us prove that there is no local minimum of macroscopic energy at fixed angular momentum and helicity. 
To that purpose, we consider perturbations of the form 6a = r(j)'j^^ and 6(^ = B'j^^(j)'j^^ /r. The corresponding 
stream function is 6(j) = (I^'mn / -^'mn- Consider first the even solutions, including the continuum solutions and the 
even eigenmodes. In that case, we have {6a) = {r(j)'MN) — 0' {i^^) = ^^'^ {4^b(I)mn) + ^{'^^mn) — ^ i'^^O — 
C B B'j^ {(j) b4>'m n) ~ since (I^'mn orthogonal to The preceding relations remain valid for the odd eigenmodes 
{m^n) provided that (M^N) ^ (m^n). Therefore, these perturbations do not change the helicity and the angular 
momentum at first order. On the other hand, for these perturbations, we have 

S^F = l- (71) 

Thus, for given B and {M^N) sufficiently large i.e. B'j^js^ > 5, we have 6'^F < 0. As a result, the critical points 
of macroscopic energy at fixed helicity and angular momentum cannot be energy minima since we can always find 
particular perturbations that decrease the energy while conserving the constraints. 
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In conclusion, the critical points of macroscopic energy at fixed helicity and angular momentum are saddle points 
since we can find perturbations making d'^F positive and perturbations making J^F negative. This analysis shows 
that all Beltrami solutions are unstable. However, saddle points may be characterized by very long lifetimes as long 
as the system does not explore dangerous perturbations that destabilize them. This motivates the numerical stability 
analysis of Sec. IIVI 

Remark: Let us consider the odd eigenmode (1,2). We have seen that it can be destabilized by a perturbation 
or by a perturbation (/>^^ at smaller scale. Let us now consider the effect of a perturbation of the form = 
1 ^ Bl(j)B^ /t and 5 a = B^r(j)B^^ where (j)B^ is the first continuum mode such that {r(j)B^) = 0. The corresponding 
stream function is S(j) = (j)B^' For this perturbation, we have {6a) = B^{r(j)B^) = 0, {(,Sa) = (^{2)^^* (0^*012) = 
and (aS^) = ^12(^^12) + ^12^* (^12^^*) — since (j)[2 orthogonal to (pB^- Therefore, this perturbation does not 
change the helicity and the angular momentum at first order. For this perturbation, we have in addition 

5^F = Bl(l-^y4>k). (72) 

This quantity is negative when < B^ corresponding to case L. This implies that the eigenmode (1,2) is also 
destabilized by the perturbation S(j) = (j)B^ which is at larger scale than the perturbations (j)^^. 

IV. NUMERICAL STABILITY ANALYSIS: ROBUSTNESS OF SADDLE POINTS 

The stability analysis performed in Sec. IIII Gl has shown that all the critical points of entropy at fixed microscopic 
energy, helicity and angular momentum are saddle points. We shall now investigate their robustness by using the 
relaxation equations derived in Paper I (for a review of relaxation equations in the context of 2D hydrodynamics, see 
[29}). These relaxation equations can serve as numerical algorithms to compute maximum entropy states or minimum 
energy states with relevant constraints. Their study is interesting in its own right since these equations constitute 
non trivial dynamical systems leading to rich bifurcations. Although these relaxation equations do not provide a 
parametrization of turbulence (we have no rigorous argument for that), they may however give an idea of the true 
dynamical evolution of the flow. In that respect, it would be interesting to compare these relaxation equations with 
Navier-Stokes simulations. This will, however, not be attempted in the present paper. 

By construction, the relaxation equations monotonically increase entropy, or decrease energy, with relevant con- 
straints. Different generic evolutions are possible: (i) they can relax towards a fully stable state (global maximum 
of entropy or global minimum of energy); (ii) they can relax towards a metastable state (local maximum of entropy 
or local minimum of energy); (iii) they do not relax towards a steady state and develop structures at smaller and 
smaller scales. In the present situation, we have seen that there are no stable and metastable states. Therefore, the 
stability analysis of Sec. IIII Gl predicts that the system should cascade towards smaller and smaller scales without 
limit (except the one fixed by the finite resolution of the simulations). This is a possible regime (see top of Fig. [TQ|) 
but this is not what is generically observed in the experiments where long-lived structures at large scales are found 
(like at the bottom of Fig. [TQ|) . Here, we explore the possibility that these long-lived structures are saddle points 
of entropy or energy with relevant constraints. These saddle points are steady states of the relaxation equations. 
Although they are unstable (strictly speaking), we argue that these saddle points can be long-lived and relatively 
robust (this idea was previously developed for 2D flows in Q). Indeed, they are unstable only for certain (dangerous) 
perturbations, but not for all perturbations. Therefore, they can be stable as long as the system does not explore 
dangerous perturbations that destabilize them. Of course, the rigorous characterization of this form of stability is 
extremely complex. In order to test this idea in a simple manner, we shall use the relaxation equations and study the 
robustness of the saddle points with respect to them. 

A. Numerical method 

1. Generalities 

Our stability analysis is based on the numerical integration of the relaxation equations 

If = -x(/3^ + M<t), (73) 

t = -^(/5^+/^^ + «)' (74) 

where D and x given functions of r and and /3, ji and a evolve in time (see below) so as to guarantee the 
conservation of the invariants. 
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In the canonical ensemble, the temperature /3 is fixed and the conserved quantities are the helicity and the angular 
momentum. The Lagrange multipliers and a{t) are computed at each time so as to guarantee the conservation 
of H and /. One may check that they are solutions of the system of algebraic equations (see Paper I) 

{DOa{t) + {{xa^) + (De)) n{t) = -/? ((x^cr) + (D^^)) (75) 

{D)a{t) + {D0m = -l3{^)- (76) 

These relaxation equations are associated with the maximization problem (fTT]) provided that, at any given time, 
the distribution of angular momentum is given by Eq. ([12]) with constant /3 (see Paper I for details). By properly 
redefining the Lagrange multipliers, they are also associated with the minimization problem (fT8|) . 

In the microcanonical ensemble, the conserved quantities are H and /. In the sequel, it will be convenient to 
fix the time dependence of /3 by imposing /^~^{t) = 2 {E — E^-^{t)) at each time. Taking into account the two other 
invariants, one may check that /3(t), a{t) and ii{t) are solution of the system of algebraic equations 

{Di)a{t) + ((xa^) + {De)) ^ + ((xV^^) + (^^)) m = 0, (77) 

{D)a{t)^{D(,)^{t)^{^)m = 0, (78) 
2{E-E''-9)/3{t) = 1. (79) 

These relaxation equations are associated with the maximization problem (p!Q|) provided that, at any given time, the 
distribution of angular momentum is given by Eq. ([12]) with (3 = f3{t) (see Paper I for details). 

In the sequel we focus on the special case D = D^r^ and x — X*^"^? where and x* are constants, that allows a 
simple numerical treatement of the relaxation equations by projection along the Beltrami eigenmodes: 



G = 



^rSn(/>n, (80) 
1 

e = ^r-^Xn(/>n, (81) 
1 

Nm 

V^ = ^rpn(/>n, (82) 



where n = (m, n),0<n<iV, 0<m<M label the modes and Nm = N x M is the number of modes. In that case, 
Eqs. ([73]) and ([71]) can be transformed into a set of Nm ODEs: 

S'n -D^ [/3Sn + /iXn + Q^(r0n)] , (83) 
X'n -X* WPn + M^n] , (84) 
Pn = ^n^^n, (85) 

where By^ is such that — A(/)n + (/>n/^^ = ^n^n- Note that the constraints couple Eqs. ([83|) - (|85]) through the parameters 
/3, a and fi. To investigate the robustness of a given stationary solution, we first perturb it with a suitable perturbation 
(see below), and then follow its dynamics thanks to the relaxation equations. Two typical time evolutions are provided 
in Fig [TOl if the solution is fragile with respect to the perturbation, it will cascade to another solution (usually the 
solution of smallest scale permitted by our resolution); if the solution is robust with respect to this perturbation, it 
will eventually return to its initial unperturbed state. To quantify the robustness of a given solution, we define a 
probabilistic stability criterion by computing the probability for the solution "to escape" from its basin of attraction. 
To that purpose, we select a threshold S and compute at each time the probability of escape 



Pesc{t) = Prob 



_(r-2(^a(0+))2) 



>S 



(86) 



using Np realizations with perturbations drawn at random at t = from a suitable ensemble (see below). This allows 
us to define "statistically fragile" solutions as those for which Pesc{t) 1 when t ^ oo, the others being referred to 
as "statistically robust". In practice, the limit t ^ oo is not accessible. We thus generalize this notion to a "finite 
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FIG. 10: Time evolution of the stream function of a solution of branch 1 of the continuum. Top: in the fragile case; Bottom: 
in the robust case. 

time", by considering the asymptotic value of pesc reached at the largest time of the simulation, tmax- In addition, the 
asymptotic value of Pesc{t) provides a mean to quantify the degree of robustness of a solution. Examples are given in 
Fig- [HI foi" a fragile and for a robust solution. As can be seen, the fragile solution is fragile whatever the threshold 5. 
However, the degree of robustness of a solution depends on the threshold S. Quite naturally, the larger the threshold, 
the more robust the solution. 




time time 

FIG. 11: Time evolution of pesc made with 200 perturbations around two Beltrami states. Top: fragile solution; Bottom: 
robust solution. Two different thresholds are used: ^ = 1.1 (continuous line) ^ = 20 (dotted line). 



Remark: although the variational problems (p!Q|) . (fTTj) and ([18]) are equivalent, and all lead to the absence of stable 
equilibrium state, the corresponding relaxation equations described previously are different. Therefore, the robustness 
of the saddle points will be different in the canonical and microcanonical settings. This can be viewed as a form of 
"ensembles inequivalence" for an out-of-equilibrium situation. 



2. Perturbations 



The stability must be investigated using perturbations that rigorously conserve the integral constraints. This puts 
some conditions regarding the shape of the possible perturbations that we can use. 

In the canonical ensemble, the integral constraints are H and /. Given an initial stationary solution (cr"*", ^'*'), the 
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perturbations (^a, 5^) must obey 

(C^^) + {^^^") + {S^Sa) = 0, 
(Sa) = 0. 

One can check that this set of constraints is satisfied by any perturbation of the form 



(87) 



Sa = er 



(89) 
(90) 
(91) 



where e is the amphtude of the perturbation, io labels an even mode while ii and 12 label two different modes different 
from io such that {cf'^(I)\^) 7^ 0. Following Eqs. ([8Q|) and ([8T]) . we have set a'^ = ^rs\(j)\ and (^"^ = ^r~'^x\(j)\. In the 
sequel, we fix the amplitude of the perturbation e through the norm = (r~^((5(j)^) by imposing 



Anry 



The modes io, ii and \2 are chosen randomly according to the following procedure: i) we draw io following a uniform 
law among the N x M even modes; ii) we draw 12 following a uniform law among the NxM or NxM-\-l modes of 
the set of allowed a"^, excluding io- This mode is therefore necessarily even for solution of continuum, and often even 
for mixed solutions; iii) we draw ii following a uniform law among the 2 x N x M even and odd modes, excluding io 
and 12' This choice allows the generation of Np random perturbations with the same amplitude Ap. 

In the microcanonical ensemble, the relaxation equations conserve in addition the energy. To satisfy this additional 
constraint, we choose the perturbations according to the same procedure as in the canonical case, and then determine 
the initial value of the temperature Po = P{t = 0) in order to guarantee the conservation of the energy [35|. As 
explained previously, this amounts to taking 

f3^'=2iE-E^-\a* + 5a,C+SO)- 

In the following, we shall group the perturbations into subclasses such that perturbations of the same class have 
the same temperature /3o or, equivalently, the same macroscopic energy E^-9{a'^ + Sa^^^"^ + S(_). Note that the initial 
temperature of the perturbation differs from the temperature of the equilibrium state which is given by 



eq 



3. Parameters 



In the sequel, we focus on the stability analysis in the case L, for the first three branches of solutions, relevant for 
comparison with experiments, see Paper III [30|. Our parameters are as follows: 

• The number of modes is Nj^ = 2 NM with = 10 (radial modes) and M = 12 (vertical modes) corresponding 
to 120 even modes and 120 odd modes. The radial and vertical lengths are = 1.2 and h = 1.4. 

• The amplitude of the perturbations is Ap = 0.05. We consider Np = 200 realizations for each given stationary 
solution. 

• The parameters and x* are both taken equal to 1. The relaxation equations are integrated using an implicit 
Heun scheme. The time step is empirically chosen proportional to j3{t = 0)~^. For j3{t = 0) = 1, the time step 
is 0.02. We have checked that this time step is small enough to guarantee the numerical conservation of / and 
H (canonical case) or I and H (microcanonical case). 



4. Numerical protocol 



For any value of A on a given branch of solutions, we proceed as follows: 
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Canonical ensemble 

i) We fix the value of the temperature {3 (it remains constant during the evolution). In the sequel, we focus on 
five arbitrary values, /3 = 1000, 100, 10, 1, 0.1, chosen so as to span a wide range. 

ii) We compute the Beltrami solution (cr^,^^) corresponding to a prescribed value of A on the given branch. 

iii) We generate Np perturbed initial conditions leaving unchanged the helicity and the angular momentum of 

iv) We evolve the perturbed initial conditions through Eqs. ([83|) - (|85]) and Eq. ()75|76p for a certain amount of time 
Microcanonical ensemble 

i) We fix the value of the energy E (it remains constant during the evolution). In the sequel, it is fixed after 
arbitrary choice of five values of the temperature, /3o = 1000, 100, 10, 1,0.1, chosen so as to span a wide range. 
Once Po has been fixed, the total energy is then fixed. It can vary from one realization to the other, but does 
not vary along the evolution. 

ii) We compute the stationary Beltrami solution (cr"^,^"^) corresponding to a prescribed value of A on the given 
branch. 

iii) We generate Np perturbed initial conditions leaving unchanged the helicity and the angular momentum of 

iv) We group together the perturbations that have the same initial temperature /3q"^ = 2{E — E"^"^ ) measuring the 
initial energy of the fluctuations (equivalently, these perturbations have the same value of macroscopic energy 

iv) We evolve the perturbed initial condition through Eqs. ([83]) - ([85]) and (j77jj79p for a certain amount of time tmax- 

B. Numerical results 

1. Robustness in the canonical ensemble 

On the three branches, we computed the value Pesd^) (computed at tmax) as a function of A for different temper- 
atures T = /3~^. The results are displayed on Fig. [121 The mixed branch (vertical dipoles) and branch 2 (reversed 
monopoles) are found very robust for high threshold, and still retain a certain degree of robustness for a small thresh- 
old, with a 40 per cent probability of escape. There is no clear dependence on the temperature. This is natural, since 
temperature can be eliminated by a suitable rescaling of time (or of coefficients D and x) and redefinition of Lagrange 
parameters. The behavior on branch 1 (direct monopoles) is more contrasted and provides a very clear transition 
around the critical value Ac ~ 0.1. For A < Ac, the probability to escape is close to 1, meaning large fragility of 
the branch. For A > Ac, the branch becomes much more robust, reaching a larger degree of robustness than the two 
other branches for high threshold, while reaching the same robustness for small threshold. The different behaviors 
are summarized on Fig. [131 

Note that the value Ac = 0.1 is somewhat arbitrary. Indeed, increasing tmax further, we observed the same 
qualitative scenario, with an increased value of Ac- We also observed that over sufficiently long time, the branch 2 
tends to become unstable, past a value of the order A = 0.25. 

2. Robustness in the microcanonical ensemble 

On the three branches, we computed the value Pesc{oo) (computed at tmax) as a function of A for classes of 
perturbations with different initial temperature Tq = (3^^. Note that the initial temperature fixes the amplitude of 
the velocity fluctuations. The results are displayed on Fig. [lH For the mixed branch (vertical dipoles) and branch 2 
(reversed monopoles), the microcanonical results do not noticeably differ from the canonical results: the two branches 
are found very robust for large threshold, and still retain a certain degree of robustness for a small threshold, with 
a 40 per cent probability of escape. There is no clear dependence on the initial temperature. There is therefore no 
ensembles inequivalence for these two branches. This is not true anymore for branch 1 (direct monopoles). Indeed, 
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FIG. 12: Canonical ensemble. Pesc{oo) as function of A on the three branches (left: branch 1, middle: branch 2, right: mixed 
branch) at different temperatures 1/13: 0.001(«), 0.01( ), 0.1(*), 1(»), and 10(»). Two different thresholds are used: ^ = 1.1 
(continuous line) (5 = 20 (dotted line). 
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FIG. 13: Robustness of the three branches in the canonical case computed at a given tmax- The lines are increasingly fat with 
increasing pesc (oo), i.e. robustness. Note that the value of Ac increases with increasing tmax. 

one still observes a transition from robustness to fragility around a critical value Ac but this quantity depends on 
the initial temperature Tq: it takes a value Ac ~ 0.1 at large initial temperatures (large velocity fluctuations) and 
then decreases to for small initial temperatures (small velocity fluctuations). The difference of robustness observed 
between the two ensembles may be seen as a kind of inequivalence of ensembles at small initial temperatures. The 
different behaviors are summarized on Fig. [151 Like in the canonical case, we checked that an increase of tmax results 
in a larger fragility of the branch 1 and 2 towards small A, at a given temperature. 

Remark: note that perturbations with small initial temperature have large macroscopic energies corresponding to 
perturbations at large scales. According to Sec. IIII Gl such perturbations are less destabilizing than perturbations at 
small scales (associated with small macroscopic energies hence large temperatures). This may explain the numerical 
results. 
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FIG. 14: Microcanonical ensemble. Pesc{oo) as function of A on the three branches (left: branch 1; middle: branch 2; right: 
mixed branch) at different initial temperatures l//3o: 0.001 (•), 0.01( ), 0.1(*), 1(»), and 10(*). Two different thresholds are 
used: (5 = 1.1 (continuous line) (5 = 20 (dotted line). 
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FIG. 15: Robustness of the three branches in the microcanonical case computed at a given tmax- Left: for a low initial 
temperature; right: for a high initial temperature. The lines are increasingly fat with increasing Pesc (oo), i.e. robustness. Note 
that the value of increases with increasing tmax • 

V. DISCUSSION 

A. Generalized ensemble inequivalence 

We have studied the thermodynamics of axisymmetric Euler-Beltrami flows and proved the coexistence of several 
equilibrium states for the same values of the control parameters. All these states are saddle points of entropy but 
they can have very long lifetime as long as the system does not spontaneously develop dangerous perturbations. We 
have numerically explored the robustness of some of these states by using relaxation equations in the canonical and 
microcanonical ensembles. The dipoles (mixed branch) and the reversed monopoles (branch 2) were found to be 
rather robust in both ensembles. Furthermore, in the microcanonical ensemble there is no dependence on the initial 
temperature on these branches. By contrast, the direct monopoles (branch 1) display a sharp transition around a 
critical value Ac. The value of Ac increases with increasing integration time. In the microcanonical ensemble, this 
value also decreases with decreasing initial temperature, resulting in a difference of robustness in the canonical and 
microcanonical ensembles. This difference may be seen as a kind of "ensembles inequivalence". This is, however, 
a very unconventional terminology since it concerns here the robustness of saddle points with respect to random 
perturbations that keep the energy or the temperature fixed, over a finite amount of time. 
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B. Bifurcation scenario 

The simulations have shown that the dipole (two-cehs solution) is relatively robust for any value of the angular 
momentum. On the other hand, the direct monopole (one-cell solution) is very fragile at low angular momentum 
but becomes robust at high angular momentum. In that case, it is even more robust than the dipole. Therefore, 
increasing the total angular momentum of the flow, one expects to observe a transition from the two-cells solution 
(antisymmetric with respect to the middle plane) to the one-cell solution (symmetric with respect to the middle 
plane). This bifurcation scenario is sketched in Fig. [161 It is reminiscent of the turbulent transition reported in the 
von Karman flow [31] in which the initial two-cells flow observed at zero global rotation suddenly bifurcates when the 
rotation is large enough. Once the bifurcation has taken place, the level of fluctuation is experimentally observed to 
decrease strongly, resulting in a decrease of the statistical temperature. In our scenario, this means that the monopole 
branch is suddenly stabilized with respect to redecrease of the total angular momentum of the flow, resulting in a 
hysteresis that has also been observed experimentally. It would therefore be interesting to investigate more closely 
the relevance of our scenario to the experimental system. This is done in the next paper [30| . 
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FIG. 16: Bifurcation scenario. Left: the system starts at A = 0, on the mixed branch, in a two-cells topology (vertical dipole); 
Middle: increasing A up to Ac (red arrow), the system follows the mixed branch. One cell grows at the expense of the other, 
resulting in a shift of the mixing layer upwards; Right: At Ac, the system bifurcates towards branch 1 (more stable), resulting 
in a one-cell topology (direct monopole). In this new state, the fluctuations are much milder (empirical fact from experiments), 
thereby allowing a stabilization of the branch towards lower A. Therefore, the monopole subsists beyond this point, even after 
a redecreasing of A, depicted by the black arrow. 



C. Richardson energy cascade 

An interesting outcome of our study lies in the fate of the solutions when they are destabilized by a dangerous 
perturbation: due to the energy minimization principle, the unstable solution tends to "cascade" towards a higher 
wavenumber solution in a way reminiscent to the Richardson energy cascade of 3D turbulence (see Fig. [10]). The 
cascade stops when the largest available wavenumber is reached, since dangerous perturbations are necessarily at 
smaller scale than the achieved state. This form of energy condensation at the smallest scale may be seen as an 
interesting counterpart (in the opposite sense) of the large scale energy condensation observed in 2D turbulence via 
the inverse energy cascade process. This is the signature of the 2D and a half nature of our system, intermediate 
between 2D and 3D turbulence. 



VI. CONCLUSION 

We have characterized the thermodynamical equilibrium states of axisymmetric Euler-Beltrami flows and proved 
the coexistence of several equilibrium states for a given value of the control parameter like in 2D turbulence [26| . We 
further showed that all states are saddle points of entropy and can, in principle, be destabilized by a perturbation with 
a larger wavenumber, resulting in a structure at the smallest available scale. This mechanism is therefore reminiscent 
of the 3D Richardson energy cascade towards smaller and smaller scales. Therefore, our system is truly intermediate 



24 



between 2D turbulence (coherent structures) and 3D turbulence (energy cascade). Through a numerical exploration 
of the robustness of the equilibrium states with respect to random perturbations using a relaxation algorithm in both 
canonical and microcanonical ensembles, we showed however that these saddle points of entropy can be very robust 
and therefore play a role in the dynamics. We evidenced differences in the robustness of the solutions in the canonical 
and microcanonical ensembles leading to a theoretical scenario of bifurcation between two different equilibria (with 
one or two cells) that resembles a recent observation of a turbulent bifurcation in a von Karman experiment |31j . 
This work was supported by European Contract WALLTURB. 



Appendix A 

We show that the first zero of F{B)^ denoted 5*, is always between the first Bi and the second B2 even eigenmode. 
To that purpose, we note that if B < B'{ then {B'^^^ - B'^) > for any (m, n) so that F{B) > 0. There is no 
discontinuity of F in the interval [0, Bi[ so that there is no zero in that interval. Consider now the interval , ^2 [. 
In that interval, F is also continuous and increasing since 

T . . > 0. (92) 

Moreover, for B {B'{y , F{B) - {(j)'-[r)'^ / {Bf - B^) -00. Similarly, F{B) - W^r)'^ / {Bf - B^) +00 when 
B (^2)"- Therefore, there exists a unique value of 5* in the range jB^^Bl^l, such that F{B^) = 0. This shows 
that the first zero of F lies in between the first two even eigenmodes. This property remains true for the successive 
values of 5*^^ and the successive even eigenmodes. 
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